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We report on systematic ab-initio investigations of Co and Cr interlayers embedded in 
Fe(001)/MgO/Fe(001) magnetic tunnel junctions, focusing on the changes of the electronic struc- 
ture and the transport properties with interlayer thickness. The results of spin-dependent ballistic 
transport calculations reveal options to specifically manipulate the tunnel magnetoresistance ratio. 
The resistance area products and the tunnel magnetoresistance ratios show a monotonous trend with 
distinct oscillations as a function of the Cr thickness. These modulations are directly addressed and 
interpreted by means of magnetic structures in the Cr films and by complex band structure effects. 
The characteristics for embedded Co interlayers are considerably influenced by interface resonances 
which are analyzed by the local electronic structure. 



PACS numbers: 72.25.Mk, 73.22.-f, 73.40.Gk, 75.47.-m 
I. INTRODUCTION 

During the last years magnetoresistive effects — in 
particular the tunnel magnetoresistance (TMR) effeclP^ 
— became increasingly imp orta nt for the fast develop- 
ing field of spintronic device^^. The first industrial ap- 
plicable TMR contacts have been built using crystalline 
MgO insulators which are epitaxially grown on as well 
as coated with iron electrode^^. Fe/MgO/Fe magnetic 
tunnel junctions (MTJs) have been extensively investi- 
gated to elucidate the mismatch between theoretically 
predictecP^ and the at least one order of magnitude 
smaller measurecP^^l TMR ratios. It turned out that 
the disparity can be attributed to differences between 
idealized (in theory) and real (in experiment) samples. 
More sophisticated theories which include imperfections, 
like interface disordei^iiJ^ or roughness effectJ^^, were 
able to close the gap between experiment and theory and 
highlight the importance of perfect interfaces. 

Although other tunnel junctions, like 
CoFeB/MgO/CoFeB MTJs with their high TMR 
ratioai^l, were put into the focus of attention, Fe/MgO/Fe 
MTJs are still intensely studied. Besides the emerging 
field of spin-torque effects^S, research is ongoing in 
search of other ways to increase the TMR ratio further. 
Instead of improving the interface quality an alternative 
means is found in the specific manipulation of the 

spin-depe ndent conductances by embedding ultrathin 
interlayer^2ill2l_ 

The insertion of a single layer-wise antiferromagnetic 
(LAFM) Cr interlayer results into even-odd oscill ations 
of the TMR ratio as a function of the Cr thicknesJ^HHl. 
In this paper we report on a first-principles study of these 
transport characteristics. We discuss the origin of these 
modulations with the apparent 2 ML- wavelength as well. 



Additionally, an analysis of the electronic transport re- 
sults for Co interlayers at both Fe/MgO interfaces is pre- 
sented. These investigations were motivated by previous 
ab-initio calculationa^Sl which predicted larger TMR ra- 
tios for MgO tunnel junctions with bcc Co(OOl) leads 
instead of Fe(OOl) electrodes. Due to the fact that Co 
grows epitaxially only up to a few monolayers on bcc sub- 
strates, a question arises whether ultrathin Co interlayers 
could be alternatively used to obtain an enhancement of 
the TMR ratios in Fe/MgO/Fe MTJs. To answer this 
question we computed the conductances and TMR ra- 
tios for small Co interlayer thicknesses and analyzed the 
results by means of the electronic structures. 

II. THEORETICAL BACKGROUND 

Our computational approach is divided into two steps. 
Firstly, the electronic structures of the MTJs are cal- 
culated from first principles. Secondly, the electronic 
transport properties are computed, using the potentials 
obtained in the first step. 

The electronic structure is determined self-consistently 
within the framework of density-functional theory 
(DFT) using a scalar-relativistic screened Kor ringa- 
Kohn-Rostoker (KKR) Green function technique^32Zl, 
The spherical site potentials were treated in the atomic 
sphere approximation (ASA) using the local spin den- 
sity approximation (LSDA) for the exchange-correlation 
potentiaP^'. Throughout this work a parameterization 
following Vosko, Wilk, and Nusaii^^ was used. 

Since structural information of Fe/MgO/Fe MTJs 
with embedded ultrathin Cr and Co spacers are not 
reported so far, we resort to a geometry of planar 
Fe(001)/MgO/Fe(001) junctions which has been de- 
termined experimentally by surface x-ray diffraction 
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analysed^fllMJ, xhis str ucture has been used in pre- 
vious theoretical studieJ^^'^. In detail, a super- 
cell geometry with six MgO layers sandwiched by 20 
Fe layers was used to compute self-consistently the 
atomic potentials. The insertion of x magnetic in- 
terlayers in planar Fe(001)/2:(Cr)/6(MgO)/Fe(001) and 
Fe(001)/a;(Co)/6(MgO)/a;(Co)/Fe(001) junctions was 
achieved by replacing x Fe monolayers (ML) at the 
Fe/MgO interfaces by x Cr or Co layers. This procedure 
implies that Cr and Co atoms occupy the same positions 
as the replaced Fe atoms; worded differently, the inter- 
layers follow the bcc structure of the Fe(OOl) leads. 

Due to the broken translational invariance in transport 
direction (z, i.e. [001]) and the in-plane translational in- 
variance, the eigenstates of the electrodes are labeled by 
in-plane wave vectors A;|| = (fc^, ky). The point group of 
the two-dimensional lattice is Amm. 

The ballistic conductance C per unit cell area A\j is 
computed for zero bias voltage in terms of transmis- 
sion probabilities (Landauer-Biittiker approach"^^) at the 
Fermi energy Ep, 
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T(fc||,£;F)dfc||. 



(1) 



The transmission probability T{k\\, Ep) is obtained by 
means of a Green function approactf^. The integration 
over the two-dimensional Brillouin zone (2BZ) requires 
typically about 90 000 fc||. The use of special fc|| point^^D 
reduces that number to 1/8. The resistance area product 
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Figure 1: (Color) (a) Dependence of spin-resolved P (ma- 
jority, minority) and AP transmissions on a logarithmic 
scale versus Cr thickness x for Bloch states at fe|| = in 
Fe(001)/a:(Cr)/6(MgO)/Fe(001) MTJs, x = 0,...,7. The 
green line is a fitted exponential to the majority transmis- 
sion. The deviation of the majority transmission from this 
fit is shown in panel (b). These data are fitted by cosine 
functions (fixed at a; = 1) with 2 ML (blue) and 2.17ML 
(red) periods, (c) Magnetic profiles of MTJs with Cr layer 
thicknesses x — 0, . . . , 7 ML. The green area highlights the 
magnetic moments of the Cr layers. 
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as normalized quantity, is used to compare the theoretical 
with experimental data. 

The optimistic TMR ratio is obtained from the RAs 
which are computed for the parallel (P) and the anti- 
parallel (AP) alignment of the two Fe(OOl) lead magne- 
tizations. 



TMR: 



RA^ 



(3) 



For the normalized TMR ratio the denominator is re- 
placed by RA^'^ + RA^. Since interfaces determine con- 
siderably the transport properties, transmittance maps 
which display T{k\\,Ep) versus fc|| need to be interpreted 
by means of the local electronic structure, rather than by 
the electronic structures of the bulk electrodes. The for- 
mer is obtained from the layer-resolved Bloch spectral 
density (SD) 
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Nai{E,k\0 = —lmTTG+{E,k\0, 



(4) 



of atom a in layer I. G'^i{E, A;||) is the site-diagonal Green 
function of that site. The trace involves integration over 
the ASA sphere and summation over spin-angular quan- 
tum numbers. 



III. RESULTS 

A. Cr interlayer in Fe(001)/MgO/Fe(001) 

In the following we present results of the thickness de- 
pendence of both the conductances and the TMR ratios 
for ultra-thin Cr interlayers which are embedded at a sin- 
gle interface in Fe(001)/x(Cr)/6(MgO)/Fe(001) MTJs. 
The Cr thickness dcr is varied in steps of monolayers 

(ML), X = 0, ...,7 with (icr,lML = 2.35 A, dcr,2ML = 

4.04 A and dcr^^^ML = rfcr,2ML + X ■ 1.69 A. 

We start with the tunneling behavior of Bloch states 
at F (fell =0). The associated transmission probabilities 
for X ~ 0, . . . ,7 are plotted in Fig. [iji. Both the minor- 
ity spin contribution for P and the AP contribution stay 
almost constant, whereas the majority spin contribution 
of P decays exponentially with an oscillatory modula- 
tion. Its decay rate is estimated by an exponential fit, 
exp(— 2«;a:) with k = 0.75 (green line in Fig. [T^). 

It turns out that an oscillatory modulation of Tj^aj 
shows up for all fc|| within the 2BZ. These oscillations 
are also present in the conductance C which is an inte- 
gral over the transmission probabilities, cq. ([T]); hence, 
there is no destructive interference which would lead to 
(complete) cancellation. Thus, it is essential to eluci- 
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Figure 2: (Color online) Fermi surface cross sections in the 
(100) plane of commensurate AFM Cr (c-afm). The nesting 
vectors along [001] (arrows) are shown in an extended zone 
scheme and are listed in Table IT] 



Table I: Nesting vectors of commensurate AFM (c-afm) Cr 
along [001], as given in Fig. [5] are characterized by their os- 
cillation periods (in ML). 

c-afm L' L" S A B C D' 15^ 
A[ML] 11.03 9.21 2.12 2.82 3.62 7.21 13.26 6.96 



date the underlying mechanism of these oscillations. To 
strengthen the discussion we focus in the following on the 
transmissions at F. 

The oscillation period can be estimated by fitting co- 
sine functions to the deviation of Tj^g^j(A;|| — 0) from the 
averaged exponential decay (black diamonds in Fig.[I|D). 
The fit with a period of 2 ML (blue) reproduces only the 
peak positions but deviates significantly in amplitude. A 
second fit, with 2.17ML period (red), hits the data best. 

The 2 ML oscillation can be explained by the local 
magnetic structure of the layer-wise antiferromagnetic 
(LAFM) Cr interlayers (Fig. The magnetic mo- 

ments of the Cr layers at the Cr/MgO interface is siz- 
ably enhanced due to the nonmagnetic MgO film. As a 
consequence of these partially uncompensated magnetic 
moments, the Cr interface layers act as spin filters for 
the tunnel currents. Hence, the latter are increased (de- 
creased) if the magnetic moments within the Cr interface 
layers are parallel (anti-parallel) to the magnetization of 
the opposite Fe electrode. Due to the LAFM growth of 
the Cr interlayer the tunnel current characteristics should 
exhibit signatures of 2 ML oscillations. The maxima of 
these oscillations should arise for P (AP) magnetic config- 
urations of the Fe leads at even (odd) multiples of the Cr 
interlayer thickness x. This behaviour was already found 
for LAFM Mn interlayers in Fe(OOl) /a;(Mn) /Vac/Fe(001) 
MTjJ^. However, the mismatch in Fig. [ija appears 
like an undersampling which cannot be satisfactorily ex- 
plained by means of the spin filter effect. The better 
match of the other oscillation period with 2.17 ML points 
to another effect which additionally influences the elec- 
tronic transport. 

There are two possible mechanisms that may explain 
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Figure 3: (Color) Complex bandstructure of bulk bcc Cr in 
the commensurate AFM (c-afm) phase along [001] for fc|| = 0. 
The real part of k± and the exponential decay rate k (imag- 
inary part of k±) are shown in panels (b) and (a, c), respec- 
tively. The color code of the bands indicates the irreducible 
representations of the point group 4mm of the associated 
Bloch states. The green circles at the Fermi energy mark the 
complex Ai-band of second kind, to k± = (0.922, 0.503)'^/"! 
which governs the transmission of majority electrons at f (cf. 
Fig.[lji). 



other wavelengths than the 2 ML. Firstly, one could 
think of spin-density waves within the Cr interlayers. 
Spin-density waves are found for Cr bulk systema^Zl and 
are related to nesting vectors of the Fermi surface. Nest- 
ing vectors that come into question are shown in Fig. [2] 
for a cross section of the Fermi surface in the (100) plane. 
The corresponding wavelengths of these vectors along the 
[001] direction (i. e. in transport and growth direction) 
are given in Table |l] The only vector that exhibits a 
wavelength which is comparable to that of the transmis- 
sion (2.17ML) is S, with an oscillation period of 2.12 ML. 
However, since we are interested in the oscillatory onset 
at fell = 0, the vector S cannot explain our findings be- 
cause it is offset from f . 

The oscillatory exponential decay of T^^^{k^^ = 0) 
is explained most promisingly in terms of the complex 
bandstructur^2ll of the Cr interlayers. Since a (continu- 
ous) dispersion relation is not defined for thin films, due 
to lack of translational invariance, we refer to the com- 
plex bandstructure of bulk Cr along [001]. The latter is 
decomposed with respect to the irreducible representa- 
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Figure 4: (Color online) Majority (P, bottom row), minority (P, middle row) and AP (top row) transmission r(fe||, Ep) within 
two-dimensional Brillouin zones for Fe(001)/a;(Cr)/6(MgO)/Fe(001) MTJs with Cr layer thicknesses of a; = 0, . . . , 7ML (from 
left to right). The two-dimensional Brillouin zones cover the range between — vr/a and vr/a. 



tions of the point group 4mm (Ai, A5, A2 and A2') of 
the associated Bloch states (Fig. [s]) . 

A complex band structure of a periodic system is the 
conventional band structure extended to Bloch vectors 
(fc|| , k±) with complex wave numbers k±. The associated 
bands can be cast into four categorieJ^ (i) real bands 
which correspond to the conventional band structure and 
have Imfc_L — 0; (ii) imaginary bands of the first kind 
have Refc_L = and Imfcj^ 7^ 0; (iii) imaginary bands of 
the second kind with Refcj^ = '^/a and Imfc_L 7^ 0; and 
(iv) complex bands with Rc/cj^ ^ 0, Re/cj^ ^ "^/a and 
Imfc_L 0. 

The imaginary part of k± is denoted as k and repre- 
sents a measure for the decay rate of evanescent Bloch 
state^3. At the Fermi energy, a complex band of the 
second kind shows up at k± — (0.922, 0.503)'^/a (circles 
in Fig.|3|. The corresponding decay rate of k = 0.78/ML 
matches well the estimated exponential decay of the ma- 
jority transmission {k = 0.75/ML, Fig. 1^). Due to the 
nonvanishing real part the exponential decay exhibits an 
oscillatory envelope with a wavelength A = '^/Rc[k±] w 
2.17 ML which agrees well with the fit in Fig. [TJ). We 
conclude therefore that the thickness dependence of the 
transmission for the majority states at the F point is very 
likely governed by this Ai state, provided the electronic 
structure of ultra-thin Cr films is well described by that 
of bulk Cr. 

Transmission maps display the spin-resolved transmis- 
sion for P and AP versus fc|| for each Cr thickness x 
(Fig. [4]). As observed for T^^^ik^ = 0) one finds within 
the majority transmission maps a clearly visible modula- 
tion of TmajC^ll) f^or fc| I -points in the center regions with 
maxima (minima) at even (odd) multiples oi x. It is 
reasonable to assume that the oscillations of these trans- 
missions are caused by the same bandstructure effect as 
it was discussed for T^aj ^ . Furthermore, the major- 
ity RA products in Fig. [5] exhibit even-odd oscillations 



as well, indicating constructive superposition of the os- 
cillations of the individual Tj^ajC'^ll)- 

Although the cloverleaf-like structures within the mi- 
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Figure 5: (Color online) Cr thickness dependence 
of (a) optimistic and (b) normalized TMR ratios in 
Fe(001)/x(Cr)/6(MgO)/Fe(001) MTJs. Since variations of 
the Cr layer thickness cannot be ruled out in experiment, 
a model with resistors in parallel connection is assumed 
to mimic Cr-thickness fluctuations (line styles indicate the 
weight w; see text), (c) Resistance area product RA for par- 
allel magnetic (P: majority, minority) and anti-parallel mag- 
netic configurations (AP), shown on a logarithmic scale. 
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Figure 6: (Color online) Transmission T{k^^, Ef) for the P (bottom row) and AP (top row) magnetic configuration within 
two-dimensional Brillouin zones for Fe(001)/3;(Co)/6(MgO)/ a;(Co)/Fe(001) MTJs, x = 0, ...,6ML (from left to the right). 
The panels on the right hand side are for Co(001)/6(MgO)/Co(001) MTJs with infinite bcc Co leads. The Brillouin zones cover 
the range between —n/a and n/a. 



nority transmission maps in Fig. [4] exhibit slight varia- 
tions which are in anti-phase to the majority transmission 
modulations, the corresponding RA products in Fig. [5] 
reveal no signatures of such an even-odd characteristic. 
The weak thickness dependence can be understood with 
help of the complex bandstructure in Fig. [3] From first- 
principles investigations on Fe/MgO/Fe it is knowrPSM] 
that A5 states are the main carrier within the minor- 
ity transport channel. Due to the A5 real band at the 
Fermi energy it is very likely that these Bloch states 
just propagate undamped through the Cr interlayer. The 
minority RA is therefore marginally affected by the Cr 
spacer and gives contributions which are similar to those 
of Fe/MgO/Fe MTJs. 

In contrast to the above finding, the pronounced modu- 
lations within the AP transmission maps in Fig. |4]— with 
maxima (minima) at odd (even) x — lead to an even-odd 
oscillation of the corresponding RA. Due to the spin- 
filter effect of the Cr interface layer RA^^ is in anti-phase 
to RAf,,^:. This behavior results in TMR ratios which 

maj 

exhibit 2 ML oscillations with periodic changes of the 
sign. The amplitudes of the TMR ratios, are with about 
7000 % and 9000 % at a; = ML and 2 ML, considerably 
larger than for the other thicknesses, with values between 
about —100% and -1-100%. This even-odd oscillation of 
the TMR ratio as a function of the Cr thickness has been 
observed experimentall5l231but with a phase shift of 1 ML. 
In more recent experiments^ it has been found that this 
phase shift depends on whether the Cr interlayer is grown 
after or before growth of the MgO barrier. However, a 
large maximum of the TMR ratio for small Cr thicknesses 
does not show up in both growth conditions. Instead, an 
exponential decay of the optimistic TMR ratio for in- 
creasing X is reportecP^. Since variations of the Cr layer 
thickness cannot be ruled out in experiment, we assume a 
model with resistors in parallel connection to mimic Cr- 
thickness fluctuations. The resistance of the mean thick- 
ness X is weighted by w, the contributions from x — 1 and 
x + 1 are weighted by (100 % — w)/2, respectively. With 
already large central weights of w = 96 % and w = 90 % 



this model is able to reproduce the principal experimental 
TMR characteristica^S. For a detailed analysis of the ef- 
fect of structural imperfections, however, one has to rely 
on more sophisticated computational approaches, like the 
coherent potential approximation or a supercell method. 

B. Fe(OOl)/ x(Co)/ 6(MgO)/ x(Co)/ Fe(OOl) 

In this section we discuss the effects of Co interlay- 
ers embedded into Fe(001)/MgO/Fe(001) MTJs on the 
electronic transport. With respect to previous theoreti- 
cal investigations of MgO barriers with bcc Co lead^^^ we 
specifically studied the thickness dependence of ultrathin 
Co interlayers which are inserted at both Fe/MgO inter- 
faces. Since Co grows epitaxially only up to few ML on 
a bcc substrate, the Co thicknesses dco are restricted to 
X < 6 ML. The monolayer separations are taken identical 
to the case of Cr interlayers. Substituting all Fe atoms 
with Co atoms the effect of semi-infinite bcc Co leads is 
studied in addition. 

First, we recall previously reported conductances and 
TMR ratios in MTJs with equal thickness variations of 
both Co interlayer^2]. The corresponding fc|| -resolved 
transmissions versus x are shown in Fig. [6] for both mag- 
netic configurations (P, AP). The shapes of the trans- 
mission maps for P configuration (bottom row) are — 
beside the case of the semi-infinite Co leads — very sim- 
ilar. In contrast to this weak dependence, one observes 
a rise and a decline of the transmission probabilities in 
the central regions of the Brillouin zones for the AP con- 
figuration for small x (top row). The maximum shows 
up at a; = 2 ML. Based on these observations it is rea- 
sonable to expect a relatively constant behavior of the P 
conductance and a maximum at x — 2 ML for the AP 
conductance. 

The elsewhere published conductances and TMR 
ratioS^ are inserted into Fig. [7] as red boxes. In accor- 
dance with the findings for bcc Co leadS^Sl^ t^^q identify 
specific Co thicknesses {x = 3 ML and 5 ML) that ex- 
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Figure 7: (Color) Co thickness dependence of (a) optimistic 
TMR ratio, P (b), and AP (c) conductances in symmetric 
Fe(001)/2;(Co)/6(MgO)/x(Co)/Fe(001) MTJs. Red squares, 
taken from Ref. |421 are for data which exhibit hotspots in 
the transmission probabilities T{k\\, Ef). The black sym- 
bols show the data with these hot spots being removed (see 
text). Panel (d) displays the layer-resolved spectral density 
(in states/Hartree) of a As-minority interface resonance with 
fc|| = (0.095, 0.008)'^/a. The blue areas mark the position of 
the 2 ML thick Co interlayers within the MTJ. 



hibit larger TMR values than those obtained for pure 
Fe/MgO/Fe junctions. In particular, the TMR ratio fol- 
lows the even-odd type characteristic vifhich shovifs up for 
the conductance C^. But, this behavior of does not 
reflect the weak thickness dependence as expected from 
the transmission probabilities (Fig. |6|. It turned out 
that this even-odd change is considerably affected by sin- 
gle interface resonances whose transmission probabilities 
contribute with up to 70 % to the conductance. These 
hot spots within the transmission maps occur preferably 
in ideal, symmetric MTJs at zero bias voltagd^i*. They 
are strongly diminished by breaking the symmetry, for 
instance by means of a tiny bias voltage or structural 
imperfections of the sample. A bias voltage of 0.02 V is 
sufficient to destroy the resonant states and to reduce 
the corresponding transmission probabilities by several 
orders of magnitude. 

The total number of hot spots in each transmission 
map of Fig. [6] is less than 10. Instead of removing the 
resonances by a small bias voltage, the transmissions of 
these states are identified and neglected. These filtered 
data are shown as black diamonds in Fig. [7] The effect 



of the resonances shows up mainly for C^. For the lat- 
ter one obtains, after an initial decrease of about 50 % 
from X = OML to IML, the expected nearly constant 
behavior. 

The filtered AP conductances agree well with unfil- 
tered ones, in particular the maximum at x — 2 ML. 
Inspecting the corresponding transmission map (Fig. |6|, 
this maximum is attributed to enhanced transmissions 
in the center region of the 2BZ which are caused by 
As-minority interface resonances within the Co interface 
layer (Fig. [7[i). In contrast to the transmission reso- 
nances (hot spots), these transmissions are unaffected by 
small bias voltages. 

Since the thickness dependence of is weakened by 
the filter procedure, that of the TMR ratio is as well. 
In particular, the maxima at x = 3 ML and 5 ML do 
not show up. The minimum of approximately 200% at 
X = 2 ML is, however, still present; it is explained by 
the increase of — up to a value comparable with 
— due to the As-minority interface resonances. For 
increasing x the TMR ratio increases monotonously up 
to a maximum at a; = 5 ML, with 4600% considerably 
smaller than the 7000% obtained for pure Fe/MgO/Fe 
MTJs. 

We note in passing that the enhanced TMR ratio 
in MTJs with semi-infinite bcc Co leads, reported in 
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Figure 8: (Color) Optimistic TMR ratios (a) and 
resistance-area products RA for P (b) and AP 
(c) of symmetric (black) and asymmetric (red) 
Fe(001)/2:(Co)/6(MgO)/t/(Co)/Fe(001)MTJs. Hot spots 
were filtered out (see e. g. Fig. [7| . 
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Ref. are not reproduced. Previous investigations of 
Fe/MgO/Fe systems have shown that conductances de- 
pend strongly on details of the calculations, in particular 
on the atomic positions in the interface region. Thus, one 
is lead to attribute the above mentioned discrepancy to 
computational details. 

Since increased TMR ratios are not obtained by sym- 
metric MTJs, we studied the effect of asymmetric MTJs, 
with one interlayer thickness fixed to 5 ML, in addition. 
The fixed thickness of 5 ML is chosen with respect to the 
largest TMR ratio in symmetric MTJs (Fig. [t]). Both 
the RA products and the TMR ratios of the asymmet- 
ric junctions do not differ significantly from those of the 
symmetric junctions (Fig. |8|. For both kinds of MTJs, 
the TMR ratio follows the characteristics of RA^^ be- 
cause RA^ is almost constant. The presence of the min- 
imum at X — 2 ML additionally substantiates that the 
As-minority interface resonances, which cause the drop 
of RA^^ , are marginally affected by symmetry breaking 
of an ideal MTJ. 



one hand a spin-filter effect which is induced by the en- 
hanced magnetic moments of the Cr interface layers and 
on the other hand the presence of complex bands which 
are formed within the Cr interlayers. The oscillations are 
therefore mixtures of 2 ML oscillations of magnetic ori- 
gin and superpositions of the individual modulations of 
the tunneling Bloch states, which can be traced back to 
the corresponding complex wave vectors. Our results fur- 
ther indicate that spin-density waves are of minor impor- 
tance for understanding of electronic transport through 
Fe/MgO/Fe MTJs with ultrathin Cr interlayers. 

The embedding of Co interlayers at both interfaces 
does not lead to an increase of the TMR ratios with re- 
spect to Fe/MgO/Fe MTJs. Please note that the refer- 
ence values for a; = (i.e. no Co interlayer) have been 
calculated under the assumption of ideal interface struc- 
tures and are therefore probably considerably overesti- 
mated. Thus, we suggest to include the effects of imper- 
fect interfaces, which are unavoidable in real samples, in 
a future investigation. 



IV. CONCLUSION 

In summary, the computed conductances and TMR ra- 
tios exhibit the oscillatory decays for thickness variations 
of ultrathin Cr interlayers, which have been found in ex- 
periment. The analysis of the associated transmission 
probabilities reveals that the tunneling of Bloch states 
is affected by the interplay of two mechanisms. On the 
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